Interrupted coarsening in a driven kinetically constrained Ising chain 
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We introduce a driven version of the ID kinetically constrained spin chain |2lj. In its original 
undriven version, this model shows anomalous coarsening following a quench to a low temperature, 
with an equilibration time that diverges as ~ exp(l/T 2 ) for T — > 0. We show that driving of constant 
rate 7 interrupts coarsening and stabilises the chain in a state analogous to that of a coarsening 
chain of age I/7. We present an analytical theory for this steady state, and demonstrate it to be in 
excellent agreement with our simulation results. 



PACS: 05.20.-y; 05.40.-a; 05.70.Ln; 64.70.Pf 
I. INTRODUCTION 

Glassy systems relax very slowly at low temperatures. 
They therefore remain out of equilibrium for long times 
and exhibit ageing the time-scale for response to an 
external perturbation (or for the decay of correlations) 
increases with the "waiting time" i w since the system was 
quenched to the low temperature, and thus eventually far 
exceeds the experimental time-scale. Time translational 
invariance is lost. 

As a result of this dynamical sluggishness, glassy sys- 
tems are highly susceptible to steady external driving, 
even when the driving rate 7 is small. (One exam- 
ple of 7 is shear rate in a rheological system.) Typi- 
cally, such driving interrupts ageing and restores a time- 
translationally invariant (steady) non-equilibrium state 
in which the time-scale defined by the inverse driving 
rate plays a role analogous to the waiting time t w of the 
ageing regime E3j. This scenario was first investigated 
in the context of neural networks 19 ll[] and has sub- 
sequently been reproduced in the diffusion of a particle 
in a random potential pfj| , ^2| ; in driven mode-coupling 
equations of the mean field p-spin mode l [q]; and in a 
driven version of Bouchaud's trap model [ p5|j29| , |3l| ] . The 
results of the latter two studies were separately exploited 
to propose a general framework for the study of "soft 
glassy materials" [p^p|J^, p^ , p6| | in which intrinsic rheo- 
logical ageing is interrupted by driving (shear straining) 
and loading (shear stress) |J. 

In this paper, we introduce a driven version of another 
glassy model: the ID Ising chain with an asymmetric 
kinetic constraint. The original undriven model (intro- 
duced by Jackie and Eisinger |2l]]) shows anomalous slow 
coarsening ( "ageing" ) following a quench to a low tem- 
perature T « 1, with an ergodic time that diverges as 
r ~ exp(const/T 2 ) for T — > 0, as solved exactly by Sol- 
lich and Evans |30| . Particularly attractive features of the 
model are i) that its glassiness emerges as a direct result 
of the dynamical constraint (without the need for any un- 
derlying assumption of quenched disorder) and ii) that it 
contains explicit spatial interactions while being simple 
enough to allow analytical progress. In what follows, our 



central result will be that steady driving interrupts coars- 
ening and stabilises the chain at an apparent age 0(1/7), 
consistently with the scenario described above. 

Ageing and driven glassy systems in general violate 
the equilibrium fluctuation dissipation theorem (FDT) 
f27f . Remarkably, however, in the ageing limit i w — > 
00 of many mean-field [|l3 12 and some non-mean-field 
HJUjH models, a non-trivial modified FDT emerges 
and defines a non-equilibrium temperature T c g [fl5|| . 
Cugliandolo et al. [fl5| proposed that an equivalent tem- 
perature should apply in the driven limit 7 — > 0, and that 
T e s(j — ► 0) and T e s (t w — > 00) should coincide. To date, 
however, the evidence H[l for this is rather limited. An 
added motivation for the present paper, therefore, is that 
the model defined here can be used (in future work (|]) 
to test this scenario further via comparison of its FDT 
predictions with those of the undriven model JlO) ] . 

The paper is structured as follows. In section || we 
summarise the results of |}0| for the coarsening dynamics 
of the undriven chain. We then define the novel driving 
rules in section III . We present simulation and analytical 
results for the steadily driven chain in sections IV and |v| 



respectively, demonstrating that driving of rate 7 halts 
coarsening at an effective age I/7. We summarise and 



give an outlook for future work in section VI 



II. COARSENING IN THE UNDRIVEN MODEL 

The model consists of a chain of L Ising spins Si E 0,1 
(1 < i < L) in a uniform field of unit magnitude, which 
is oriented such that a spin's energy Ei(si) — s^. Peri- 
odic boundary conditions apply: the left neighbor of si 
is sl. The dynamics are subject to the following con- 
straint: at any time, only those spins whose left neigh- 
bor is up are allowed to flip. For these "mobile" spins, 
the rate of down-flips is 1 while the rate of up-flips is 
e — exp(— 1/T). The equilibrium distribution is unaf- 
fected by this constraint: detailed balance is obeyed and 
the static distribution is the trivial one prescribed by the 
Hamiltonian H — X)f=i s i- ^ n contrast, the dynamics 
are rendered very slow at low temperatures, for which 
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the equilibrium concentration of up spins (that facilitate 
the dynamics) c = e/(l + e) is small. 

In this section we review [[30|, which solved exactly 
the relaxation of a chain prepared at a low temperature 
Tel via rapid quench at time t w = from a high initial 
temperature Tj = oo. At any time i w during the relax- 
ation, the authors of j3(]] described the system's state 
using the concept of "domains" . As shown by the verti- 
cal lines in 

- -- 1 1 - - 1 1 

(in which the spins are represented by "•" for clarity) a 
domain is defined as an up-spin, and all the down-spins 
separating it from the next up-spin to the left. 

Immediately after the quench, the average domain 
length d = 1/c = 0(1). The system then relaxes to- 
wards the low temperature's equilibrium state in which 
d = 0(l/e), and thus in which, in the limit e — > 0, there is 
zero probability of finding an up spin in any finite length 
of chain. Hence, in this limit (e — > at fixed chain length) 
the down-flipping of spins is irreversible, and the relax- 
ation comprises a coarsening process in which adjacent 
domains progressively coalesce with one another. 

The mechanism for this coalescence is as follows. Con- 
sider, at some time after quench, a domain of length 
d <C 1/e together with its left-bounding up-spin. For 
the purposes of the present argument, we assume that 
the left-bounding spin is "clamped" and consider how 
the right-bounding up-spin relaxes. Because of the con- 
straint, before the relaxation can occur a facilitating up- 
spin has to be generated immediately to the left of this 
spin, via a propagation of the up state rightwards from 
the left-bounding up-spin. The relaxation is thus im- 
peded by an energy barrier, the height of which is the 
maximal number of spins that are ever up within the 
(original) domain at any instant during this relaxation 
process. A central result of p(J] is that for domain lengths 
2 n ~ 1 < d < 2 n this barrier scales as n, leading to a re- 
laxation time-scale 0(e~ n ). 

Hence in the limit e — > the dynamics comprise well 
separated stages, the n th of which has time-scale e~" and 
results in domains of index n being destroyed by coales- 
cence with their right neighbors. In logarithmic time, 
v w = — log/j w /loge = Tlog(i w ), the n th stage collapses 
to the point v w = n. The average domain size d thus 
exhibits step-wise increases at successive integer values 
of v w , as seen in figure 1 of |3C|| . 

Within this coarsening regime (d <C 1/e), the full do- 
main length distribution P(d) can be calculated using 
an exact independent interval treatment which states 
that no correlations can build up in the length of ad- 
jacent domains provided none are present in the initial 
state. At stage n the distribution obeys |30| 

d T P(d,r) = ]T P(d-d',T)[-d T P(d',r)] (1) 

2"- 1 <d'<2" 



in which the rescaled time r = te n can take any positive 
value r > in the limit e^O. Equation |l| describes the 
coalescence of the "active" domains of length d' < 2™ 
with neighboring domains of length d — d'. Its initial 
condition is the domain length distribution at the end of 
stage n — l of the dynamics, denoted P n (d) = P(d,T—>0). 
Using generating functions it can be shown that 

d-l 

P n+ i(d) = P n (d) +Y,Pn(d~ d')P™\d') 

d'=l 

1 d-l 

--J2 p n ct (d-d')p: ct (d>) + ... (2) 

d' = l 

in which P n act (d) is the active part of the distribution 
(zero for d > 2") and in which... denotes a series of 
convolutions of increasing order. Equation |^ holds only 
for d > 2": all active domains (d < 2") disappear in the 
n th stage. 

The weight of the distribution shifts to larger d at each 
stage of coarsening. See figure 2 of |50|. A scaling limit 
is approached for large stage number n: the re-scaled 
distribution P n (x — d/2 n ~ 1 ) = 2 n ~ 1 P n (d) converges to 
the limit P agc (a;) which obeys, for x > 2, the scaling 
counterpart of equation ||: 

ipage (|) = Page(z) + £ dx' P agc (x - X>)P£(x>) 

-~ f dx'Pg(x - x')P^{x') + ... (3) 

In |30]| the exact solution for P agc (x) was shown (via the 
re-summed Laplace transform of |3|) to be 

p ^) = 2l—^ 1 n— 6 (£*.-* ■ 

m— 1 r— 1 \s — 1 / 

(4) 



III. DEFINITION OF THE DRIVEN MODEL 

In this section we incorporate non-Hamiltonian driving 
into the model. As a preliminary step, though, we rede- 
fine the relaxational dynamics slightly, extending state 
space such that each spin Si S —1,0,1 and re-defining 
the Hamiltonian H = X)i=i \ s i\- The uniform field has 
thus been replaced by a potential well for each spin, with 
a minimum at = 0. As before the dynamics are con- 
strained: only those spins for which the left neighbor has 
a value 1 or —1 are allowed to flip via the usual thermal 
processes in which the transition rate for Sj : 1 — > and 
for Sj : — 1 — ► is 1, and for Sj : — > 1 and Sj : — > — 1 
is e. 

So far, of course, the model can be exactly mapped 
onto the original one by a trivial relabelling Sj = — 1 — > 
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Si = 1 and re-scaling e — > e/2. Our motivation for in- 
troducing the — 1 state is to make a loose analogy with 
glassy rheological models [ ^TUTt|JT^ | in which a local state 
of high energy (here |sj| = 1) can have either positive or 
negative local stress (here \st\ = +1, |sj| = —1). If we 
define a global stress a — \ Ylt= l s «> the — 1 state allows 
a state of macro scopically zero stress, which still has in- 
ternal local stresses (some positive, some negative) and 
a non-zero rate of internal dynamical rearrangements. 

We now incorporate steady driving into this three state 
version. Loosely this mimics, in a stochastic way, the 
standard rheological experiment of applying shear strain 
of constant rate j. To do this, we impose a flip rate of 7 
for Si : — 1 — * and 7 for Sj : — ► 1. This is additional to 
the constrained rates defined above and free of the kinetic 
constraint. The driving rates ui for this extended model 
can therefore be summarised as follows: 

u(si : 0— v 1) = |sj_i|e + 7 
u(si : — > -1) = |si_i|e 

u(si : 1 -> 0) = |si_i| 
u>(si : -1 ->• 0) = |si_i| + 7, 

(5) 

in which periodic boundary conditions impose so = s l- 

This stochastic straining clearly tends to increase the 
global stress, as required intuitively. We note, though, 
that our stochastic rules only make sense for 7 > 0. For 
negative 7 we would redefine the driven contribution to 
the rates as equal to 7I for the transition s,; : 1 — > and 
|7| for Si : — » —1. The model is in this sense singular 
at 7 = 0. 



IV. SIMULATION RESULTS FOR THE 
STEADILY SHEARED MODEL 

We simulated the driven chain using a waiting time 
Monte Carlo technique combined with a binary search 
algorithm for locating the mobile spins, following [ p0[ . 
For each run we initialised the chain either in equilibrium 
(with 7 = 0) at a low temperature T = — 1/loge <C 1, 
or by quenching to T from T — 00. For the quenched 
case, we then let the system relax according to the un- 
driven rules (7 = 0) until a start-up time t s , when we 
set 7 to the non-zero, constant value of interest. Equi- 
librium initialisation formally corresponds to a quenched 
chain subsequently allowed to relax until t s — 00, and 
in this case we applied the non-zero 7 from the start of 
the simulation. In order to explore the hypothesis that 
driving restores a steady state analogous to the state of a 
coarsening chain of age t w = 1/7, we chose values of I/7 
corresponding to the waiting times studied in |30fl for the 
undriven chain. Specifically, therefore, we are interested 
in the low temperature limit in which e — > 0, 7 — > 0, at 



fixed values of v = log(7)/ log(e) = Tlog(l/7) that are 
large compared to 1 (weak driving), but small enough 
that the system remains far from equilibrium. (To avoid 
possible confusion we note that v becomes large as 7 be- 
comes small, since 7 = e v with e< 1. Although this is 
at first sight a counter intuitive way to characterise the 
shear rate, we chose this particular definition for v as 
the closest possible analogy that of v w for the undriven 
chain.) 

In each run we monitored the stress a and total energy 
E as functions of time t w . Results for the quenched ini- 
tial condition with start-up time t$ = are shown in fig- 
ure U for e = 0.01 and various values of v. As v becomes 
larger (7 smaller) we can make the following observa- 
tions. At early times, E and a have time evolutions that 
are independent of 7, and that (we have checked) are the 
same as those of an undriven chain. In contrast, after a 
crossover at time 0(1/7), E and a approach steady-state 
values. We have checked by repeating the simulation for 
the different initial conditions described above that these 
steady values do not depend upon the initial state or (for 
the quenched case) the start-up time t s . For the remain- 
der of the paper, we shall be concerned only with the 
ultimate steady state, and not the kinetics of its forma- 
tion. 

The steady state stress is re-plotted in figure [2] as a 
function of — v (which increases with 7) for various small 
values of e <C 1. (In rheological parlance, (7(7) is the flow 
curve.) It appears to be approaching a step-like func- 
tion as temperature is tracked towards zero, comprising 
plateaux separated by jump-wise discontinuities at inte- 
ger values of v. 




FIG. 1. Stress a (solid lines) and energy E (dashed lines) 
vs. scaled time v w — — log(t w )/ log(e) following a quench at 
time tw = 0, with driving commenced also at t s — 0; for all 
curves e = 0.01. The parameter v = log 7/ log e has values 0.5 
(O), 1.0 (□) 1.5 (O) 2.0 (A) and 2.5 (v)- Each curve was 
obtained from a single run for a chain of length L — 2 1B . 

For all the steady states studied, we found the concen- 
tration of —1 spins to be small (O(e)), and hence that (to 
within such corrections) the stress, the energy, and the 
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FIG. 2. Steady state flow curves plotted on a log scale vs. 
— v = — log 7/ log e. Simulation results for three values of 
e = exp(l/T) are shown, obtained from a single run for a spin 
chain of length L = 2 16 . Bold line: theoretical prediction for 
T ^ 0. 



concentration of +1 spins coincide (consistently with the 
results of figure |l|). Using the usual domain description 
(in which for definiteness we neglect the —1 spins, taking 
only 1 spins to constitute the domain boundaries), we 
have then a = 1/d where d is the average domain length. 

We note the striking similarity between the depen- 
dence of logcZ upon scaled waiting time v w = Tlog£ w 
in the coarsening chain (figure 1 of J30|), and the de- 
pendence of log a = — log d upon the scaled driving rate 
— v = — Tlog(l/7) in the steadily driven chain (figure || 
of this paper). This already gives us a strong indication 
that the steady state of a chain driven at rate 7 is anal- 
ogous to the state of a coarsening chain of age t w = 1/7. 

This hypothesis is confirmed by our simulation data 
for the full domain length distribution P(d) on the de- 
veloping plateaux of the flow curves. In particular, we 
find that the distribution for the steadily driven chain 
at a given v is closely analogous to that of a coarsening 
chain for is w ~ v. For a given v (or v w ), both display 
a discontinuity at the same v (z/ w ) dependent cutoff and 
have very similar averages. Both shift abruptly to larger 
values of d as v (^ w ) crosses successive integers, but are 
unchanged as v is swept between integers. In or- 
der to maintain the closest notational analogy with |30| , 
we denote by P n (d) the distribution P(d) in the limit of 
small e for values of v such that n — 1 < v < n. Our 
results for P n {d) for v — 0.5, 1.5, 2.5, corresponding to 
n = I, 2, 3, are shown in figure || and are, as just noted, 
very similar to the counterpart results of figure 2 of |3C| l 
for the coarsening chain. 

In the coarsening chain, the discontinuous shift of P(d) 
as is w crosses successive integers arose from the waiting 
time crossing the time-scales for successive coarsening 
stages. Likewise in the driven chain it arises from the 



inverse driving rate crossing these same time-scales. (See 
the theory section |v| for more details.) 

These simulation results therefore show that (consis- 
tently with the phenomenology of other driven glassy 
models) driving interrupts coarsening and stabilises the 
chain in a state that is strikingly analogous to that of a 
coarsening chain of apparent age 0(1/7). 
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FIG. 3. Domain length distributions P„(d) on the plateau 
v = n — 1/2 of the flow-curves. Open symbols and lines: theo- 
retical results for n = 1 (squares), 2 (diamonds), 3 (triangles). 
Full symbols: simulation results for a chain of length L = 2 16 
and e = 10~ 4 (n = 1,2) and e = 10 -3 (n = 3). Inset: scaled 
predictions 2 n_1 P„(d = 2 n ~ 1 x) vs. x for n = 1, 16. 



V. THEORY 

We now present a theory aimed at calculating the 
steady-state distribution P n (d) just defined. We assume 
from the outset that no correlations exist between the 
lengths of adjacent domains. We also assume (consis- 
tently with our above remarks) that the concentration 
of —1 spins is negligible. We will return to justify these 
assumptions in more detail at the end of this section. 

Consider, then, a chain with an O(l) population of 
+1 spins and a complementary O(l) population of 
spins. In the absence of driving, the only process affecting 
the chain would be the coarsening described above, with 
down-flipping of 1 spins occurring irreversibly (for fixed 
chain length and e — > 0), and the average domain length 
increasing as domains of length 2" _1 < d < 2™ are de- 
stroyed on a time-scale 0{e~ n ) by coalescence with their 
right neighbors. 

In the steadily driven chain, this coarsening process 
is balanced by intra-domain driven up-flipping of spins 
Si : — > 1. For the finite domain lengths to which we 
shall restrict ourselves, this driven domain intersection 
occurs on a time-scale 0(1/7) = 0(e _ly ). For a non- 
integer value of v such that n — 1 < v < n, this time- 
scale sits between the time-scales e~(™ _1 ) and e~" for the 
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adjacent coarsening stages n— 1 and n (recall section |J), 
and separates from them in the limit e — > 0. 

On this driving time-scale, therefore, any domain of 
index n' < n — 1 (present either as a remnant of the ini- 
tial condition, or as a result of a driven up-flip — * 1 
within an existing domain of length d > 2™ _1 a distance 
d < 2 n ~ l from the domain boundary) must relax in- 
finitely quickly. We thus expect lim^o Pn{d) = for 
such domain lengths d < 2 n ~ 1 = d c . This defines the 
cut-off length d c observed in the simulation data above. 

Domains labelled by n' > n — 1, i.e. of length d > d c , 
on the other hand, coarsen infinitely slowly on the driving 
time-scale: the only dynamical processes that can affect 
these "long" domains are those initiated by driven up- 
flip of an intra-domain spin. This up-flip can occur at 
a distance that is cither < d c or > d c from either end, 
giving 4 separate cases: 

• A. If the original long domain was longer than 2d c + 
1, and if the up- flip occurred at a distance of at 
least d c + 1 from both ends, we see creation of two 
shorter domains that are still both "long" in the 
sense that d > d c : 



■1- 



(For definiteness all the diagrams in this section 
assume an arbitrary cutoff value d c = 4.) 

Considering the chain as a whole, such processes 
lead to the destruction of domains of length d > 
2d c + 1 at a rate 



j(d-2d c -l)P n (d) 



(6) 



and to the creation of domains of length d > d c at 
a rate 



27 ]T P n (d') 
d'=d+d c + l 



(7) 



B. If the up-flip occurred within a distance d c of 
the left-hand end of the original domain, but at a 
distance greater than d c + 1 from the right hand 
end we see a process such as 

...1 1 1... 

I 

...l-l 1 1... 

I 

...1 1 1... 

which does not need to be considered further since 
it results in no net change. 



• C. If the up- flip occurred within a distance d c of 
the right-hand end of the original domain, but at a 
distance greater than d c from the left hand end, the 
"short" (right) sub-domain will relax immediately 
by coalescence with its right neighbor: 

...1 1 1... 

I 

...1 1... 

I 

...1 1 1... 

leaving two "long" (d > d c ) domains, with the 
boundary between them shifted to the left. 

Processes such as this lead to the creation of do- 
mains of length d at a rate 



d-\-d c d c 

7 £ P n(d ')+7£ 

d'=d+l d' = l 



--d c +d' + l 



P n (d-d') 
(8) 



and to the destruction of domains of length d at a 
rate 

7 9(d - 2d c )d c P n {d) + 7 6(2d c - d+ l){d - d c - l)Pn(d) 

2d c oc 

+ 7 P„(rf) E (d 1 -d c ~l)P n {d') + d c J2 P ^ d ') 

_d'=d a + l d'=2d c + l 



where the discrete Theta function is defined by 

0(rt — to) = 1 for n > to 
= for n < to. 



(9) 



• D. If the original domain was of length d < 2d c , 
both the sub-domains could be short (d < d c ). For 
this class as a whole (i.e. averaging over the po- 
sition of the up-flipped spin) there are two sub- 
classes of possible outcome, which by symmetry 
each occurs with a probability one half. First, the 
freshly flipped up-spin could relax before the right- 
bounding spin of the original domain and we would 
see no net change. Alternatively, we could see 

...1 1 1... 

I 

...1 ■ -1 ■ -1 1... 

I 

...1 • -1 1... 

I 

...1 1... 

which is essentially "aided coalescence". For the 
chain as a whole, processes such as these lead to 
destruction of domains of length d at a rate 
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jP n {d)^(2d c -d + l) 
+ 2 P n(d) A f Pn(d')(2d c -d' + 1) 

and creation at a rate 

| ^ (2d c -d' + l)P n (d')P n (d-d'). (10) 

(i'=d c + l 

Combining all these processes we get an evolution equa- 
tion dtP n {d) = ... for the "long" domains d > d c , which 
we set equal to zero (steady state) and solve numeri- 
cally using an iterative procedure. The solutions for 
n = 1, 2, 3 for which d c = 1, 2, 4 respectively are marked 
as open symbols figure [| and give excellent agreement 
with the simulation results. We also used these solu- 
tions to calculate the stress, a = 1/d. As expected, this 
exhibits discontinuous jumps at integer values of v as 
the driving time-scale crosses successive coarsening time- 
scales and P(d) shifts ^continuously to larger d. It is 
marked as the solid line in figure |^ and again agrees ex- 
cellently with the simulation data. 

As d c — > oo, a scaling limit d c P{d) = P(x — d/d c ) 
is approached. See the inset of figure ||; and figure ||. 
Taking the limit d c — > oo at fixed x = d/d c in the steady 
state equation just derived, we find that this scaling state 
must obey the equation 

= -g(x)P(x) + 2 dx'P{x') - / dx'P(x') 

J X J X 

+ f dx'f(x')P(x-x'). (11) 
Jo 

We note for use below that the function f{x) (thus de- 
fined) is discontinuous at i = 1, and the first deriva- 
tive g^ x \x) of g is likewise discontinuous at x = 2. 
In principle, equation |ll| contains all the information 
needed to calculate P analytically. For the counterpart 
state P a ge (a:) in the undriven chain, the closed expression 
(equation [|) was found ||(| as the self-consistent solution 
of a simple algebraic relation between the Laplace trans- 
forms G and H of P a ge(^) and P agc (x)Q(2 — x) respec- 
tively. The corresponding transform of equation |ll| for 
the driven case is a complicated differential relation be- 
tween G, H, G (1) and and we have been unable find 
a self-consistent analytic solution. However our numeri- 
cal results (see figure ^) demonstrate that P(x) is (as ex- 
pected) very similar to its undriven counterpart P& ge (x): 
both have a unit Heaviside discontinuity at x = 1 and 
show similar decay for x > 1. Although the discontinuity 
in the first derivative Pigl(x) at x = 2 (strongly appar- 
ent in the dashed curve of figure |4|) is less noticeable in 
the solid curve for the driven state P, it is revealed by 
numerical differentiation in the inset of figure m 



It was shown in |50| that the ageing scaling distribu- 
tion P agc has a finite discontinuity in its k th derivative 
at x = k + 1 for all integer k > 0. We have already seen 
numerically that the driven scaling distribution shares 
the discontinuities for k = 0, 1. We shall now outline an 
analytical argument which can be used to show that in 
fact the driven distribution shares all of these discontinu- 
ities. We confine ourselves to k > 1 since our analysis has 
already captured the discontinuity in P itself at x = 1. 
Differentiating equation [n] once we get 

- g{x)P w {x) = g (1) {x)P{x) + P{x) + P(l + a;) 

- f dx' f{x')P {1 \x-x'). (12) 
Jo 

On the right hand side (RHS) of this expression, g^ x \x) 
is discontinuous at x = 2 (as noted above), while the 
integral over the infinite discontinuity §(x — x' — 1) of 
the differential P^'(x — x') in the last term's integrand 
gives 0(3 — x)f(x ~ 1), which is discontinuous at x = 
2. All other terms on the RHS are continuous for all 
x > 1. Hence P^(x) has a finite discontinuity at x = 2. 
By performing successive differentiations, we can extend 
this argument to arbitrarily high k: p( k \x) has a finite 
discontinuity at x — k + 1 for all k > 0. 




ti i i i i i i i i i i i i i r^i 

0.5 1 1.5 2 2.5 3 3.5 4 



_ x 

FIG. 4. Scaling distribution P(x) for the driven chain (solid 
line) and P a _ sa (x) for the coarsening chain (dashed line). In- 
set: derivative dP/dx for the driven distribution, showing a 
slight discontinuity (too small to be discernible in the undif- 
ferentiated data of the main figure) at x = 2. 

Our central result can therefore be summarised as fol- 
lows. We have shown that a steadily driven chain ap- 
proaches a scaling state in the limit of small 7. We have 
shown that this scaling state is strikingly analogous to 
the counterpart scaling state of a coarsening chain of age 
t w = I/7. In particular, both states have a unit Heavi- 
side discontinuity at x = 1 , and (more generally) a finite 
discontinuity in P^(x) at x = k + 1. 

Despite the striking qualitative similarities in the dis- 
tributions P(x) and Pa g e(aO, there are obvious quanti- 
tative discrepancies. This is not inconsistent with the 
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phenomenology of other glassy models: the ageing state 
of Bouchaud's trap model (25), for example, is analogous 
but not identical to its steadily driven counterpart |3l| . 
It is nonetheless instructive to consider the origin of the 
discrepancies. An obvious candidate is our introduction 
in the driven model of the —1 spin state. However, in 
steady state the population of such spins is small. Indeed, 
we have checked that domain distributions produced by 
driven simulations without the —1 state agree with those 
that include the —1 state to within 0(e). 

The discrepancies must therefore be of dynamical ori- 
gin. During coarsening, the only process is domain co- 
alescence. Therefore information can only propagate up 
the distribution (to larger d) . In this way the scaling limit 
Page(x) is completely determined by P(x') for x' < x, as 
seen in equation]^. In contrast, in the driven chain P(x) 
is connected to all x' — > oo, since driven domain inter- 
section acts on domains of all lengths; see equation [n]. 
Furthermore, in the ageing chain the domain lengths set 
the dynamical rates only via the coarsening time-scales 
e~™. In contrast, the driving dynamics are further sensi- 
tive to domain length through the additional feature that 
a longer domain is more likely to be intersected. 

The most apparent difference between P age (x) and 
P{x) is the much smaller discontinuity in PW(x) at 
x = 2. This can be explained as follows. The discon- 
tinuity in P^(x) arises from an interplay between the 
discontinuity in g^(x) at x — 2 and in f(x) at x = 1 
(recall equation |l2]) . Physically, g(x) is the rate at which 
domains of scaling length x are destroyed. This is discon- 
tinuous because driven intersection of a domain of scaled 
length x < 2 can produce two active domains (scaling 
length < 1), whereas only one can be active for x > 2. 
Similarly, the discontinuous part of f(x') results from the 
production of domains of length x via aided coalescence 
(process D above) which can only occur for x < 2. These 
two effects are of comparable (not identical) magnitude 
but opposite sign resulting in only a small discontinuity in 
P^(x). In contrast, the discontinuity in Pagl(x) arises 
from a single process, the origin of which can be seen 
from equation ||. The second term on the RHS describes 
the production of domains of length x via coalescence of 
an active domain of length x' with a right neighbor of 
length x — x' . However if x — x' is itself active, the right 
neighbor can relax before the x' domain, resulting in a 
domain of length > x. The rate of producing domains 
of length x is thus reduced, as encoded in third term on 
the RHS. The cutoff in the active distribution at x = 2 
means that the derivative of this term is discontinuous 
at x = 4. When transferred to the next coarsening stage 
(encoded in the LHS of the equation), this discontinuity 
appears in Pigo(2). The higher convolutions not shown 
in equation ^ are smooth enough not to affect this argu- 
ment. 

We finally return to justify our two assumptions: 



firstly, that the population of —1 spins is negligible in 
the steady state. As noted above, we have already nu- 
merically observed this population to be 0(e) (small); we 
are now in a position to show this theoretically, as fol- 
lows. Over a time-scale ~ 1/7, the only mechanism in 
which —1 spins can be created involves propagation of 
the facilitating state (—1 or 1), via constrained thermal 
activation, to within a distance d < d c to the right of 
one of the existing facilitating spins. Denoting by r c (d) 
the time-scale on which domains of length d coarsen, the 
time-scale upon which a —1 spin is created at a distance 
d from an existing facilitator is r c (d)/e. Once created, 
such a spin will relax back to the state on a time-scale 
T c (d). Hence in the steady state, the population of — 1 
states must be O(e), which is indeed small in the limit 
considered here. 

Our second assumption was that of "independent in- 
tervals". This was used in |3^] for the coarsening pro- 
cess of the undriven chain, and in that case is prov- 
ably exact. We have not been able to prove its strict 
validity for this driven case. However, our simulation 
results show that the relevant correlation function C = 
E"=o(^ " (d))(d 3+ i (d»/E?=o(4f ~ (d)) 2 (m which 
dj is the length of j th domain from the left hand end of 
the chain, and n is the total number of domains) is not 
greater than 10 -3 in any of the steady states considered; 
the assumption therefore is likely to be reasonable. Note 
that our numerical observation of a non-zero correlation 
function could still be consistent with the independent 
interval approximation being exact in the limit e — » at 
fixed v. Indeed, we have observed that, as e is tracked 
downwards, the numerical value of the correlation func- 
tion gets smaller. 

VI. SUMMARY AND OUTLOOK 

We have incorporated driving dynamics into the kinet- 
ically constrained spin chain of [pl[| , and presented sim- 
ulation results showing that the coarsening dynamics of 
the undriven chain (as solved in f30|) are interrupted by 
steady driving. Consistent with the broader glassy liter- 
ature |l|,|l^,^l]j20|j3|J|,||,|l| we have found that steady 
driving stabilises the chain at an apparent age set by the 
inverse driving rate. We have presented a theory for this 
steady state, demonstrating it to be in excellent agree- 
ment with our simulation results. We have shown that 
a scaling state is approached at small 7, and that this 
scaling state has very similar properties to its counter- 
part scaling state reached at long times in the coarsening 
chain. 

We now outlook some possible directions for future 
work. Above, we focused on a chain which is steadily 
driven at a constant rate 7. We also noted the loose 
analogy of this driving scenario to that of constant shear 



7 



rate in a rheological system. In the spirit of this rheologi- 
cal connection, we can identify possible analogues of two 
other standard rheological tests - step strain and step 
stress - which it would be interesting to investigate fur- 
ther. For a step strain of size 70 we promote, at the time 
i w of strain application, a fraction 70 of —1 spins (chosen 
randomly) according to Sj : — 1 — > and of spins (again 
chosen randomly) Si : — ► 1. We do this without regard 
to the kinetic constraint. For all other times the system 
merely relaxes under its undriven constrained dynamics. 
(Note that the step at t w is just the "impulsive limit" 
of the above steady shear case: jdt = 70, with dt — > 
and 7 — > 00.) For a step stress of size 00, we apply the 
same dynamics just defined for the step strain up until 
the time i+. (We can merely rename 70 by ctq because 
the "spring constant" k = 1.) For t > t w we implement 
the "constant" strain-rate dynamics defined above, but 
with 7 continuously adjusted to ensure that <7o remains 
(on average) a constant. 

As noted in the introduction, it would also be interest- 
ing to study FDT in the driven steady state (at constant 
7) to see if any effective temperature emerges, and (if so) 
whether it coincides with any effective FDT temperature 
of a coarsening chain of age t w = 1 /j. This is the subject 
of a forthcoming publication Q . 
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